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ABSTRACT 

In the framework of turbulence dynamo, flow motions amplify a weak seed 
magnetic field through the stretching of field lines. Although the amplification 
process has been a topic of active research, less attention has been paid to the 
length scales of magnetic field. In this paper, we described a numerical study on 
characteristic lengths of magnetic field in magnetohydrodynamic turbulence. We 
considered the case of very weak or zero mean magnetic field, which is applicable 
to the turbulence in the intergalactic space. Our findings are as follows. (1) At 
saturation, the peak of magnetic field spectrum occurs at ~ L Q /2, where L is the 
energy injection scale, while the most energy containing scale is ~ L Q /5. The peak 
scale of spectrum of projected, two-dimensional field is ~ L . (2) During the stage 
of magnetic field amplification, the energy equipartition scale shows a power-law 
increase of ~ t 1 ' 5 , while the integral and curvature scales show a linear increase. 
The equipartition, integral, and curvature scales saturate at ~ L , ~ 0.3L , and 
~ 0.1 5Lo, respectively. (3) The coherence length of magnetic field defined in the 
Faraday rotation measure (RM) due to the intergalactic magnetic field (IGMF) is 
related to the integral scale. We presented a formula that expresses the standard 
deviation of RM, ctrm, in terms of the integral scale and rms strength of the 
IGMF, and estimated that o~rm would be ~ 100 and ~ a few rad m~ 2 for clusters 
and filaments, respectively. 
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Introduction 



It is well established that the universe is permeated with magne tic fields (see, e.g., 
Kronberg|l994j ). and yet the origin of them is not well understood (see iKulsrud fe Zweibel 



2008 



and references therein). The problem of cosmic magnetism can be divided into two 
parts - the origin of seed fields and their amplification. In this paper, we are concerned with 
the latter. 

If a weak seed magnetic field is introduced into a turbulent medium, flow motions stretch 



field lines and ampli 



see, e^ 



Batchelor 



V the field. Such turbulent ampl i fication has been studied since 1950's 
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Rvu et al.ll2008l : 



Cho et al. 



20091 ). The amplification goes through three stages: (1) Eddy motions in turbulence are 
fastest on the smallest scale, which is the dissipation scale. Therefore, the stretching of 
magnetic field lines occurs most actively at the dissipation scale first, and the magnetic 
energy grows exponentially. (2) The exponential growth stage ends, when the magnetic en- 
ergy becomes comparable to the kinetic energy at the dissipation scale. The follow-up stage 
is characterized by a linear growth of magnetic energy, and a gradual shift of the peak of 
magnetic field spectrum to larger scales. (3) The amplification of magnetic field stops, when 
the total magnetic energy becomes comparable to the kinetic energy. A final, statistically 
steady, saturation stage is reached. 



Ryu et al.l (120081 ) proposed a scenario in which turbulent flow motions are induced via 



the cascade of the vorticity generated at cosmological shocks during the formation of the large 
scale structure. Based on a model of turbulent amplification, they estimated the strength 
of the intergalactic magnetic field (IGMF); (B) ~ a few fiG in clusters, ~ 0.1/xG around 
clusters, and ~ 10 nG in filaments. In their model, the intracluster medium is at a stage 
close to the saturation one. But the IGM in filaments are still at the linear growth stage. 



The IGMF has be en observe d with Faraday rotation measure (RM) ; for instance, I Clarke et al. 



(j200lh in clusters and Ku et all (120061 ) outside clusters. For the IGMF with (B) = 0, the 
mean value of RM is expected to vanish, and the standard deviation represents the observa- 



tion (see, e.g.. lRyu et al.lll998l ). Hence 



crm — 0.81 n e B\\ rms I 



L 7 radm- 



(1) 



where n e , B\\, and / are in units of cm -3 , /iG, and pc, respectively, has been used to extract 
the strength of the IGMF. Here, B\\ rms is the root-mean-square (rms) strength of line-of-sight 
magnetic field, I is the coherence length, and L is the path length. But the formula has been 
applied without a clear definition of the coherence length. 
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Here, we study length scales of magnetic field in magnetohydrodynamic (MHD) turbu- 
lence with very weak or zero mean magnetic field, intended for application to the IGMF. 
Several characteristic lengths, which have been introduced in previous works, are considered. 
Characteristic lengths of initial seed fields in the early universe could be either very large or 
very small, depending on how the seed fields were generated. However, the field structure at 
the linear growth and satu ration stages is i ndependent of the seed fields as long as they are 
sufficiently weak (see, e.g.. lCho et al.ll2009l ). Therefore, the characteristic lengths should be 
determine d by the energy injection scale (or the outer scale) only. And based on a scaling 
argument, ISchekochihin fc Cowleyl (120071 ) stated that in the linear growth stage, the scale 
at which the stretching is most active grows as l s oc t 3 ^ 2 . 

In this paper, through numerical simulations, we investigate the growth and saturation 
of characteristic lengths. Then, we show how the coherence length in the RM due to the 
IGMF is defined in ter ms of the charact eristic lengths, and present a formula for <Jrm- And 
based on the model of iRyu et al.l (120081 ). we estimate urm for clusters and filaments. 



2. Simulations 



We solved the incompressible MHD equations in a periodic box of size 2ir using a 
pseudo-spectral code. We drove turbulence in Fourier space. The forcing function has the 
following form: f(x,t) = Y^j=x ffoj) exp(zfcj • x) + complex conjugate, where f(kj) is a 
complex vector that is perpendicular to the wavevector kj. The 22 forcing components are 
nearly isotropically distributed in the range 2 < k < Vl2, where k = \k\. The phase of 
each forcing component randomly fluctuates, but it has a correlation time of approximately 
unity. The amplitude of each forcing component also randomly fluctuates. On average, 
each forcing component injects similar amount of energy B In physical space this forcing 
f(x, t) corresponds to statistically homogeneous driving on large scales. Only the solenoidal 
component of the velocity field was driven^. The strength of energy injection was tuned, 
so without magnetic field the rms velocity becomes unity, V rms ~ 1, at saturation. In this 
representation, V can be viewed as the velocity measured in units of the rms velocity of the 



2 The wavenumbers used are 2 (3 components), (12 components), 3 (3 components), and vl2 (4 
components). Therefore the peak of energy injection occurs at ko y 2.5, and the en ergy injection scale is 
Lq ~ 2.5. A similar discrete sampling of forcing can be found in iTao et al.l (|1993l ). where they used 13 
components. 

3 Since we deal with incompressible turbulence, we applied solenoidal forcing, 
cases, the properties of turbulence may depend on the nature of forcing (see, e.g., 
for compressible turbulence). 
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system. The density is unity and the magnetic field is multiplied by y An in our simulations, 
so B can be viewed as the Alfven speed in the same units. The magnetic field consists of 
the uniform background field and the fluctuating field: B = B$ + b. At t = 0, the magnetic 
field had either weak uniform component (when Bo 7^ 0) or only random components (when 
Bq = 0; Run 256H8-B 0, see below), and the velocity had a support between 2 < k < 4 in 
the wavevector space. We conside red only the cas e wher e the kinetic viscosity , u, is equal 
to the magnetic diffusivity, r\. See ICho Vishniad (120001 ) and ICho et al.l (120091 ) for further 
details of simulations. 

Simulations are denoted with XY-B Z, where X = 256 or 512 refers to the number 
of grid points in each spatial direction, Y = H8 or P refers to hyper (and their order) or 
physical viscosity/diffusivity, and Z refers to the strength of mean magnetic field. In this 
paper, Runs 256H8-B 0, 256H8-B 10~ 3 , 512P-B 10~ 3 , and 512#8- J B 10~ 4 are discussed. 
In 256H8-B 0, the mean field was zero, and initially the spectrum of magnetic field peaked 
at k ~ 70. All the runs have either very weak or zero mean magnetic field, since we intend 
to apply the results to turbulence in the IGM0, as noted in Introduction. Besides the initial 
magnetic field configuration, Runs have either different numerical resolution or different 
viscosity/diffusivity to explore their effects. In 512P-i?olO _3 physical viscosity/diffusivity 
was used, while in others hyperviscosity/diffusion was used to extend the i nertia l rangqj. 
Runs with 256 3 grid points are a subset of the simulations used in lCho et al.l (120091 ). For the 
pr operties of t urbul ence with different initial magnetic fields, please refer an extensive study 



m 



Cho et al.l (120091 ). Figure 1 shows the time evolution of V 2 and B 2 in four simulations. 
Here, the kinetic energy and magnetic energy densities are V 2 /2 and B 2 /2. Although the 
four simulations have different set-ups, their time evolution looks similar. We can clearly 
see three stages of magnetic energy evolution. In this limit of week or zero mean field, the 
resulting turbulence is globally almost isotropic. 



4 The origin of seed magnetic fields in the IGM can be either cosmological or astrophysical. Theories 
in favor of cosmological origin suggest that weak seed magnetic fields were created in the early universe or 
during the structure formation era. Our initial magnetic field intends to mimic the cosmological origin of 
seed fields. 

5 It is well known that hyper-dissipation causes the bottleneck effect, unphysical flattening of energy 
spectrum near the dissipation range. However, since the bottleneck effect is negligible at small wavenumbers, 
we believe the effect does not alter the shape of the magnetic energy spectrum at small wavenumbers. Most 
length scales discussed in this paper rely on the shape of magnetic energy spectrum at small wavenumbers. 
Therefore, we believe the use of hyper-dissipation does not affect our results much. Indeed, our results show 
that the length scales do not strongly depend on numerical resolution or forms of dissipation. 
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3. Various Length Scales of Magnetic Field 

Peak scale of spectrum of magnetic field, Le{h)'- Figure 2a shows the spectra of 
velocity and magnetic field at an epoch of saturation. In E v (k), the peak occurs at the 
energy injection scale. Note that the peak appears at k = 2 rather that at 2.5, which is a 
numerical artifact of discrete binning of k. In Ej,(k), the peak occurs at a smaller scale. 

Largest energy containing scale of magnetic field, LkE(k)'- With E(k) oc 1/k 
representing a power spectrum that contains equal amount of energy in each decade of k, 
the peak scale of kE(k) defines the largest energy containing scale. Figure 2a shows that for 
the magnetic field, the peak of kE b (k) OCClirS clt £L SCcllc even smaller than that of Eb(k). 

Peak scale of spectrum of projected magnetic field, LiE(k)/k' When a three- 
dimensional (3D) scalar quantity is projected onto a two-dimensional (2D) plane, the energy 
spect rum of the projected quantity becomes E 2D {k) oc E 3D (k)/k (see, e.g.. lCho fc Lazarian 



20021 ). It is a bit more complicated for vector quantities, but the qualitative behavior should 
be similar. Figure 2a shows that for the magnetic field, the peak of Eb(k)/k occurs at a scale 
close to that of E v (k). Note that the largest energy containing scale of projected magnetic 
field corresponds to the peak scale of kE 2D (k), which is the same as the peak scale of E b (k). 

The time evolution and saturation of L E (k), L kE (k), and L E (k)/k, normalized with the 
energy injection scale Lq, are shown in Figure 2b. The rugged profiles are again a consequence 
of discrete binning of k. Two points are noteworthy. (1) At saturation, Leoa ~ Lq/2, 
LkE(k) ~ Lq/5, and L E (k)/k ~ Lq. (2) The growth pattern of those scales in the linear growth 
stage is not clear owing to the ruggedness in the profiles, but it seems to be between ~ t and 

Energy equipartition scale, L eq : The energy equipartition wavenumber, k eq , is de- 
fined by 



/ 

Jkr 



E v (k) dk = E b (k) dk. (2) 

Jo 



The time evolution and saturation of L eq are shown in Figure 3a. We expect that the 
stretching of magnetic field lines is most active at this scale. If th en, L pq should represent 



In (se e Introduction) and follow ~ £ 5 in the linear growth stage (ISchekochihin fc Cowley 



20071 ). Indeed, L eq grows with a power-law index consistent with the theoretical prediction 
in Figure 3a. At saturation, L eq becomes close to Lq. 

Integral scale, L int : The length scale of magnetic field (and velocity too) is often 
characterized with the integral scale, which is defined by 

jE b {k)/k dk 

Unt = 271 jE b (k)dk (3) 
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It is well known that the integral scale has the same order of magnitude as the longitudinal 
and transversal integral scales, L\ and L t , respectively, and in incompressible isotropic tur- 
bulence (with reflection invariance on average) they are related by L/ = 2L t = (3/8)L int (see 



Monin fc Yaglomlll975f ). 



Curvature scale, L curv : We also consider a typical radius of curvature of field lines, 
L C urv We define it as the distance r at which the average correlation drops to 1/e, 

(B(x)-B(x + r)) x I 
(B(x) ■ B(x)) x e, {) 
where the two points at separation r are located along the same magnetic field line and the 
average is taken over x. The factor 1/e is an arbitrary choice. 

The time evolution and saturation of Li nt and L curv are shown in Figure 3b. We cal- 
culated L curv only for runs with 256 3 grid points. Unlike L eq , L int and L curv seem to grow 
linearly in the linear growth stage (although the reason of the linear growth is not clear). At 
saturation, L int ~ 0.3L and L 

curv ~ O.lSX/o- Hence, for instance, we model the growth and 

saturation of L int as 

f (0.3/45) L x t/t eddy , if t/t eddy < 45 
int ~ \ 0.3 L , if t/t eddy > 45 _ U 

where the eddy turnover time is defined with the vorticity around the energy injection scale 
at saturation as t eddy = ^-j ^injection- 



4. Faraday Rotation Measure 

With (RM) = for the IGMF, the standard deviation of RM is 

,l N 2\ V2 

2 




crm = 0.81 n e ( I / B\\ds I ) rad m . (6) 



o 



Note that in this work, the density is assumed to be constant. Here, J B\\ds is the projected, 
2D magnetic field. 

Without loss of generality, we take x as the line of sight direction. The projected field 
can be written as § 

p2tt i-2tv p2ir N / 2 

I B\\ds= I B x dx& dx J2 B x (k)e lX - k 

&x fey fez = -~ N / 2 



6 We considered turbulence in a cubic box of size 27r, as noted in Section 2. That is, we assumed 
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N/2 



27T ^ 4(0,^,^)6^+^), 
k y ,k z =-N/2 



(7) 



where we used f Q n dx e lxkx = 2ir5o t k x - Then, the square-average of the projected field becomes 

1 



2tt - - M 

Buds 



o 



(27T) 5 



2?T /•27T 



/ 



2tt 



B x dx 



N/2 



(2tt) 2 ^ 1^(^ = 0,^,^) 

k y ,k z =-N/2 



2tt) 2 / dfc,,dfc 



(2vr) 2 / gft(fe) 



dife. 



In the last step, (\B x \ 2 )k x =o = (l/2)(\B\ 2 )k x =o statistically in the k x = plane and Ej,(k) 
Airk 2 (\B\ 2 / 2) were used, and dk y dk z was substituted with 2nkdk. Finally, using {B 2 )/2 
j E b (k)dk, 



2tt n - m 

B\\ds 



o 



(B 2 ) (2tt) 2 jE b (k)/k dk 
4 J E b (fc) dfc 

(g 2 )L mt (27T) 

4 



(9) 



So far, we have assumed that the box size is 2n. When the box size is L (or, when the path 
length is L)@, 27r's in Equation flU]) should be replaced by L. Hence, the standard deviation 
of RM becomes 



r, ni - -Brms y/ Li n tL _o 

crm = 0.81 n e rad m . 



(10) 



Since B\\ rms = B rms j V3, the coherence length used in Equation (CQ) should be given as 
I = (3/4)L int . 

Equation (TTOT) can be app lied to estimate RMs due to the IGMF. Here, we employ 
the model of iRyu et al.l (120081 ) for turbulence and magnetic field in the IGM. In clusters, 



a flow of period 2ir in the three directions of space (see iLesieurl I2008L for further details), 
case, we have a usual Fourier series expansion: B x (x) ~ Y^ulX k =~n/2 B x {k)e lX k and 



points in each side. Here, k = (k x , k y , k z ), and Xi mn 



Jinn n j x,,x,K. lmn ,^ ' k , where N is the number of grid 

denotes the coordinate of a grid point (l,m,n) in the 
computational grid. 

7 The periodic simulation bo x size can be used as a proxy for the physical path length of the system (see, 
e.g., discussion in Lesieur 2008h . 
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turbulence is near the saturation stage with t/ teddy ~ 30, where t = t age , the age of the 
universe, and t e dd y = ^rms- From Equation (j5J), then, L int /L ~ 0.2. If we take the energy- 
injection scale L ~ 100 kpc, which is approximately the scale height of cluster core, L int ~ 20 
kpc. With n e ~ 10~ 3 cm -3 , -B rms ~ a few /xG, and the path length L ~ 1 Mpc, we get 
~ 100 rad m~ 2 for clusters, which agrees with the observed RMs in clusters ( jClarke et al. 
20011 ). In filaments, on the other hand, with t/t e dd y ~ 10, turbulence is expected to be still 



in the linear growth stage, and L int /L ~ 1/15. We may take the energy injection scale 
L ~ 5 Mpc, which is the typical thickness of filaments. It is also the typical radius of 
curvature of cosmological shocks in filaments (IRyu et al.ll2003l ). which would be the major 
sources to drive turbulence there. The power spectrum of v cur i, the curl component of flow 
motions which satisfies the relati on V x v r . nr i = V x v, in the large scale structure of the 



universe peaks around 5 Mpc too (IRvu &: Kang 



120081). 

field strength in filaments quoted in lRyu et al.l (120081 ) is (B) ~ 10 nG. But we note that the 



value depends on how it is averaged. With the data of IRyu et al.l (120081 ). (B 2 ) 1 ^ 2 ~ a few 
x 10 nG, (pB)/(p) ~ 0.1 pG, and ((pB) 2 ) 1 / 2 / {p 2 ) 1/2 ~ a few x 0.1 pG, in the warm-hot 
intergalactic medium (WHIM) with T = 10 5 — 10 7 K which mostly composes filaments. The 
average value of {{pB) 2 ) 1 / 2 / (p 2 ) 1 ^ 2 should be most relevant to RM. Then, grm for filaments, 
if they intersect the line of sight with right angle, would be 



<?rm ~ 1-5 



n,, 



10 



-5 cm 3 



B r 



J int 



0.3 pGJ \300kpc 5Mpc 



1/2 



rad m 2 . 



(11) 



Normally filaments would not intersect the line of sight with right angle. Smaller angles 
result in larger path lengths and so larger orm-, and then the typical value of <jrm for 
filaments could be a few rad m -2 . We note that the values of IRMI toward the Hercules 



and Perseus-Pisces supercl usters reported i n lXu et al.l (120061 ) is an order of magnitude larger 
than the above value, and Ku et al.l (120061 ) quoted the path length, L, which is about two 
orders of magnitude larger. 



5. Summary 

Our findings are summarized: 

(1) We studied different characteristic scales of magnetic field in MHD turbulence with very 
weak or zero mean magnetic field. They saturate at ~ 0.1 — ILq, where Lq is the energy 
injection scale. 

(2) During the linear growth stage of magnetic energy, the energy equipartition scale follows 
L eg oc t 3 / 2 , while the integral scale follows L int oc t. 

(3) The integral scale (actually (3/4)Li nt ) is the relevant scale for RM. We obtained a new 
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formula for the standard deviation of RM, orm (see Eq. [10]). 

(4) We estimated a^u for filaments as well as for clusters of galaxies. 

Finally, we note that our findings are based on a small number of incompressible numer- 
ical simulations. They appear to be rather insensitive to numerical resolution and forms of 
dissipation. However, when different initial conditions or different types of forcing are used, 
some of the results may be different, a question which should be investigated in the future. 

We thank the anonymous referee for useful suggestions. JC was supported in part 
by Korea Research Foundation (KRF-2006-331-C00136). DR was supported in part by 
Korea Science and Engineering Foundation (R01-2007-000-20196-0). JC and DR were also 
supported by Korea Foundation for International Cooperation of Science and Technology 
(K20702020016-07E0200-01610). 

REFERENCES 

Batchelor, G. 1950, Proc. R. Soc. Lond. A, 201, 405 
Brandenburg, A., & Subramanian, K. 2005, Phys. Reports, 417, 1 
Cho, J., & Lazarian, A. 2002, ApJ, 575, L63 
Cho, J., & Vishniac, E. 2000, ApJ, 538, 217 

Cho, J., Vishniac, E., Beresnyak, A., Lazarian, A., & Ryu, D. 2009, ApJ, 693, 1449 

Clarke, T. E., Kronberg, P. P., & Bdhringer, H. 2001, ApJ, 547, Llll 

Federrath, C, Klessen, R. S., & Schmidt, W. 2008, ApJ, 688, L79 

Kulsrud, R., & Anderson, S. 1992, ApJ, 396, 606 

Kulsrud, R., Cen, R., Ostriker, J., & Ryu, D. 1997, ApJ, 480, 481 

Kulsrud, R., & Zweibel, E. 2008, Rep. Prog. Phys., 71, 046901 

Kronberg, P. P. 1994, Rep. Prog. Phys. 57, 325 

Lesieur, M. 2008, Turbulence in Fluids, 4th ed. (New York: Springer) 

Monin, A. S., & Yaglom, A. M. 1975, Statistical Fluid Mechanics, (Cambridge: MIT Press) 



-10- 



Ryu, D., & Kang, H. 2008, in Numerical Modelings of Space Plasam Flows: Astronum 2007, 
ASP Conf. Series Vol. 385, ed. N. V. Pogorelov et al. (San Francisco: APS) p44 

Ryu, D., Kang, H., & Biermann, P. L. 1998, A&A, 335, 19 

Ryu, D., Kang, H., Clio, J., & Das, S. 2008, Science, 320, 909 

Ryu, D., Kang, H., Hallman, E., & Jones, T. W. 2003, ApJ, 593, 599 

Schekochihin, A. & Cowley, S. 2007, in Magnetohydrodynamics: Historical Evolution and 
Trends, ed. S. Molokov et al. (Berlin: Springer), p85 ( jastro-ph/0 507686) 

Tao, L., Cattaneo, F., & Vainshtein, S. I. 1993, NATO Advanced Study Institute: Solar and 
Planetary Dynamos, p303 

Xu, Y., Kronberg, P. P., Habib, S., Dufton, Q. W. 2006, ApJ, 637, 19 



This preprint was prepared with the AAS IATgX macros v5.0. 



- 11 - 




Fig. 1. — Time evolution of V 2 and B 2 . Time is given both in units of the eddy turnover 
time defined with the vo rticity around th e energy injection scale at saturation, t e dd y = 
1/ 00 i n j ec tion, (bottom) as in Ryu et al. ( 2008h and in units of the energy injection scale divided 
by VV 2 + B 2 at saturation (top) as in bho et all (120091 ). 
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Fig. 2. — Left panel: Spectra of velocity, E v (k), and magnetic field, E b (k), at t /teddy = 59. 
The dissipation spectrum is also shown. Although hyper-dissipation affects the spectra 
near the dissipation scale, it does not strongly influence the shapes of the spectra at small 
wavenumbers. Right panel: Time evolution of peak scales of Eb(k), kEb(k), and Eb(k)/k. 
A's are the scales normalized with L Q . 
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Fig. 3. — Time evolution of the energy equipartition scale (left panel) and the integral and 
curvature scales (right panel). A's are the scales normalized with L . 



